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A Consistent Picture of a Collapsing Bose-Einstein Condensate 
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We present a consistent picture of the recent JILA experiments on the collapse of a Bose- 
Einstein condensate (BEC) [E. A. Donley, et al., Nature 412, 295 (2001)] based on a generalized 
Gross-Pitaevskii equation that incorporates the effect of a three-body loss. The questions we 
address in this paper are: (1) Why does a BEC survive, albeit partially, after the collapse? (2) 
Why does the number of the remnant BEC atoms observed at JILA exceed the critical value? 
(3) What is the origin of the atom burst? (4) What is the origin of the atom jet and is it 
coherent? We review our predictions on pattern formation in the course of collapse and on the 
split instability of vortices that can test our collapsing theory. 
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1. Why attractive interaction? 

One of the unique features of gaseous Bose-Einstein 
condensates (BECs) is their unprecedented controllabil- 
ity when applied to the creation, manipulation, or prob- 
ing of the BEC system, and as a result, many previously 
gedanken experiments have become or are expected to 
become a reality. A prime example is the dynamics of a 
collapsing BEC, which mimic some of those of collapsing 
and exploding stars. *) 

In a gaseous BEC system, the interaction between 
atoms is characterized by the s-wave scattering length a. 
At zero magnetic field, 7 Li, 2 ) 85 Rb, 3) and 133 Cs 4) have 
negative scattering lengths, which implies that atoms in- 
teract attractively and the condensates tend to collapse 
upon themselves. We may thus expect the density in- 
stability and the ensuing collapsing dynamics of the con- 
densate which are similar to those of supernova explo- 
sions. The attractive interaction also drastically changes 
the nature of superfluidity and that of topological exci- 
tations. 

In the present paper we present a consistent picture 5-7 ) 
of the so-called Bosenova experiments recently performed 
at JILA. 8 ) Our theory is based on a generalized Gross- 
Pitaevskii (CP) equation that incorporates the effect of a 
three-body recombination loss. ^Theoretical analyses on 
the Bosenova experiments using similar approaches 10-12 ) 
-as well as analyses based on elastic condensate colli- 
sions 13 ) and molecular formation 14-16 )- have also been 
reported. We shall review some new phenomena on pat- 
tern formation in the course of collapse 5 ) and on the 
split instability of vortices 17 ) that our theory predicts 
and hence can test our collapsing theory. 

2. Controlled collapse of the condensate 

It has been well established that an attractive BEC 
may exist in a confined system when the number of BEC 
atoms is below a certain critical value. 2 ' 18 ' 19 ) In this 
case, the zero-point energy due to confinement serves as 
a kinetic barrier against collapse, allowing formation of 
a metastable BEC. It is through this dynamically gener- 
ated barrier that the system is believed to tunnel when 



the number of BEC atoms is just below the critical value 
and the system collapses spontaneously. 20,21 ) 

The subject of the present paper is the controlled col- 
lapse of the condensate. 3, 8 ' 19,22 ) Recently, the JILA 
group utilized the Feshbach resonance to switch the sign 
of interactions from repulsive to attractive, and caused 
a BEC to collapse. 3 ' 8,19 ) After the switch, the atomic 
cloud shrank and eventually disappeared because the 
cloud became too small. That is, the BEC collapsed. 
Surprisingly, however, just after the collapse they ob- 
served a burst of atoms emanating from the remnant 
BEC, 8 ) a phenomenon similar to supernova explosion. 
Before this experiment was performed, there had been 
controversy as to what happens to a BEC once it be- 
gins to collapse. Some researchers predicted a complete 
collapse, while others 5,9 ) predicted a partial one. The 
JILA experiments clearly demonstrated that the collapse 
is partial. 

To understand the underlying physics of this partial 
collapse, let us consider a BEC confined in an isotropic 
harmonic trap. The dynamics of a BEC in a dilute-gas 
regime are well described by the GP equation 
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V 2 ^ + -(x 2 +y 2 + \z 2 )?P + g\iP\ 2 TP, (1) 



where ip represents the order parameter or the "wave 
function" of the condensate, normalized as J \ip\ 2 dr = 1, 
and A is the asymmetry paramter of the trap. In Eq. (1), 
the time and the length are measured in units of the in- 
verse trap frequency ui~ 1 and the natural length scale of 
the harmonic trap do = (S/mu) 1 ' 2 , respectively, where 
m is the atomic mass. The last term in Eq. (1) describes 
the mean-field interaction between atoms whose strength 
g is given in terms of the number of BEC atoms No and 
the s-wave scattering length a as g = AwNoa/do. 

Figure 1 shows the time development of the collapsing 
condensate wave function 6 ) obtained by numerically inte- 
grating Eq. (1) using a finite-difference method with the 
Crank- Nicholson scheme. 18 ) The superimposed dashed 
curves represent Gaussian fits having the same width as 
the numerical ones. At t — 0, we increase the strength 
of attractive interaction to just above its critical value 
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Fig. 1. The solid curves show the profiles of the wave functions 
\ip(r, t)\ at t = 0,2.5,2.95, and 2.96, obtained by numerically 
integrating the Gross-Pitacvskii equation (1) with A = 1. The 
dashed curves indicate Gaussian fits that have the same widths 
as the exact solutions. 



and cause the system to collapse [Fig. 1(a)]. At t = 2.5, 
the peak height grows significantly, but the wave func- 
tion fairly well maintains the Gaussian form [Fig. 1 (b)] . 
However, at t = 2.95, the peak height grows very rapidly 
and the Gaussian approximation clearly breaks down 
[Fig. 1(c)]. Remarkably, just a little after this (t — 2.96), 
we see that the peak height has grown sufficiently that 
inelastic collisions would be substantial, while in most 
other regions inelastic collisions are negligible [Fig. 1(d)]. 
This inclastic-collision-dominated regime is extremely lo- 
calized and occupies only one millionth of the condensate 
volume. Thus the condensate has an extremely small 
atom drain -or a tiny black hole- at its center. 

That the collapse occurs only in an extremely localized 
region and nowhere else in the BEC is the primary rea- 
son that the collapse occurs only partially and we have a 
remnant BEC. 

3. What value should we take for K 3 ? 

To follow the dynamics of the system after the collapse, 
we include in the GP equation the atomic loss due to the 
three-body recombination loss: 9 ) 



:^ 3 |VlV, (2) 



where K 3 is the three-body recombination loss rate co- 
efficient. (The two-body dipolar loss has little effect on 
the collapsing dynamics.) 

The collapsing dynamics depend crucially on the value 
of K 3 -the three-body recombination rate coefficient. 
However, no direct measurement of K 3 has been re- 
ported near the Feshbach resonance peak. It is there- 
fore important to devise a method for determining K 3 
from other experimentally measurable quantities. For 
this purpose, we note that the peak density in the course 
of collapse ceases to grow when the interaction term 
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Fig. 2. The mean burst energies feB^burst plotted against 
a 2olla P se/^ 3 for various values of a collapsc and K 3 , where k B 
is the Boltzmann constant and a co n apse is the value of the s- 
wave scattering length a when the collapse occurs. The dashed 
line is a theoretical fit (3) obtained in Ref. 6 ' 



is counterbalanced by the three-body loss term, i.e., 
g\rb\ 2 ~ -ftr 3 |V>| 4 /2, and hence we have \rb\ 2 - 2g/K 3 . As 
explained in the next section, the atom burst is caused by 
a slight imbalance between the kinetic energy and the at- 
tractive interaction. We may thus expect that the mean 
burst energy is proportional to the interaction energy at 
the time when those two balance, i.e., g\ip\ 2 ~ 2g 2 /K 3 . 

To check this hypothesis, we have performed exten- 
sive numerical simulations by varying g and K 3 and cal- 
culating the mean burst energy for each set of g and 
K 3 . We calculate the energy distribution of the burst at 
t = implosion + tt/2 when the burst atoms spread max- 
imally and hence can easily be distinguished from the 
remnant BEC. Figure 2 summarizes our numerical sim- 
ulations. 

From Fig. 2 we find that when T bu rst < 200nK the 
mean burst energy is indeed proportional to g 2 /K 3 oc 
a coiia P sc/^3 and is given by 



mean burst energy ~ 2.6ti 



a 

3 ^collapse 

m 2 K 3 ' 



(3) 



By taking the experimentally measured mean burst en- 
ergy ~100nK, we obtain K 3 ~ 2 x 1CT 28 cm 6 /s, which 
we will use in the following discussions. 

All relevant parameters have thus been determined ex- 
perimentally and our theory includes no adjustable pa- 
rameters. 

4. Intermittent implosions 

The solid curve in Fig. 3 shows the time evolution 
of the peak height of the wave function. 5 ) In the ini- 
tial stage of the collapse, the atomic density grows very 
slowly. However, a rapid implosion breaks out at t ~ 
2.96; a blowup of this implosion is shown in the inset. It 
is interesting to observe that implosions occur not once 
but several times intermittently, and at each implosion 
only several tens of atoms are lost from the condensate. 
This is a consequence of the fact that the collapse occurs 
in an extremely localized region. 
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Fig. 3. Time evolution of the peak height \ip(r = 0,t)\ of the col- 
lapsing wave function according to Eq. (2) with A = 1. The 
dashed curves (referring to the right axis) show the fraction of 
the atoms remaining in the trap, which is defined as the ratio of 
the number of remnant BEC atoms N IernnSin t (t) plus the num- 
ber of burst atoms TVburst (t) at time t against the initial number 
of BEC atoms iVo(0). The inset shows an enlarged view of the 
intermittent implosions. 



The physics of the intermittent implosions may be un- 
derstood as follows. 

• Initially, the attractive interaction is made to dom- 
inate the kinetic energy, so the condensate shrinks 
towards the central region and the peak height 
grows until the growth is stopped by inelastic col- 
lisions. Only the localized high density region col- 
lapsed and disappeared due to the atomic loss. 

• After the atoms at the density peak are suddenly 
removed due to the three-body loss, the remain- 
ing BEC still peaks steeply and the zero-point ki- 
netic pressure arising from Heisenberg's uncertainty 
principle barely surpasses the attractive interaction. 
Thus the central portion of the condensate subse- 
quently expands due to this extra zero-point pres- 
sure. This slight surplus of kinetic energy is the ori- 
gin of the atom burst and explains why the energy 
of the burst atoms is as low as lOOnK. 

• As soon as the atom burst occurs, the peak height 
decreases very rapidly, which makes the attractive 
interaction again dominate the zero-point pressure, 
thereby inducing the subsequent implosions. 

• This implosion-burst cycle repeats, giving rise to 
intermittent implosions, until the number of the re- 
maining atoms becomes so small that the attractive 
force no longer overcomes the kinetic pressure. 

For an isotropic trap, the number of remnant BEC 
atoms is always found to be smaller than the critical 
number iV cr ~ 0.57c?o/|a|, 18 - ) as shown in Fig. 4(a). How- 
ever, for an anisotropic trap this is not necessarily the 
case, because the burst atom cloud cannot refocus at the 
center of the BEC as efficiently as for an isotropic trap. 
For the cigar-shaped geometry used in Ref., 8 - 1 we find 
that the fraction of the remnant BEC atoms does not 
show a significant dependence on the initial number of 
BEC atoms N as shown in Fig. 4(b), in agreement with 
the observation in Rcf. 8 ' 
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Fig. 4. Dependences on the initial number of atoms JVo of the 
remnant iVremnant , burst iVburst , and missing atoms -/V m i ss i ng for 
an isotropic trap (a) and an anisotropic trap (b). [Reproduced 
from Fig. 4 of H. Saito and M. Ucda, Phys. Rev. A 65, 033624 
(2002).] 



It should be noted that the time scale of the inter- 
mittent implosions is determined by the competition be- 
tween the loss and self-focusing of the atoms, and is much 
shorter than the time scale of the inverse trap frequency 
The oscillatory behavior in the peak density as- 
sociated with intermittent implosions should be distin- 
guished from other oscillatory behaviors, such as the 
collapse-and-growth cycle 9,23 ) (which is the oscillation 
in the number of BEC atoms) and small collapses found 
in Ref., 9 ) both of which time scales are much longer than 
that of the intermittent implosions. 

Perhaps the most convincing evidence of the intermit- 
tent implosions is the gradual decay in the number of 
remnant BEC atoms (Fig. 1 of Ref. 8 )). The dotted curve 
in Fig. 5 shows the decay in the number of remnant BEC 
atoms, which is in excellent agreement with the experi- 
mental observation. 8 ) Our simulations thus reveal that 
the gradual decay is caused by a series of intermittent 
implosions. 

As shown in Fig. 5, iV rcmnant remains constant for some 
time, which is denoted as t C oiiapso, after the switch of 
interactions. Open circles in Fig. 6 show i co iiapse obtained 
by our numerical simulations as a function of 

^collapse; 

where "error bars" allow for experimental uncertainties 
±2ao, where ao is the Bohr radius, of the initial value 
Oinit = of the s-wave scattering length (see Ref. 8 ' for 
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Fig. 7. The integrated column densities seen from the radial di- 
t [ms] rection under the same conditions as the JILA experiments. > 

[Reproduced from Fig. 6 of H. Saito and M. Ueda, Phys. Rev. 
Fig. 5. Dependences on time of the number of remnant BEC A 65, 033624 (2002).] 

atoms (dotted curve) and of the number of the remnant plus 
burst atoms (dashed curve). [Reproduced from Fig. 1 of H. 
Saito and M. Ueda, Phys. Rev. A 65, 033624 (2002).] 
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Fig. 6. Dependence of t co ii apso on a co n apsc . Open circles show the 
results of our simulations, where the "error bars" allow for exper- 
imental uncertainties of the initial value of the s-wave scattering 
length, a; n ; t . 8 ) Filled circles show the experimental results of 
the JILA group. 24 ) 



details); the upper and lower end points of each error bar 
correspond to ai n j t = — 2ao and di n it = 2ao, respectively. 
Filled circles show the latest experimental results of the 
JILA group. 24 -* The agreement looks excellent. 

5. Atom Jets 

In the JILA experiments, prolonged atomic clouds, re- 
ferred to in Ref. 8 ) as "jets", were observed when the 
collapse was interrupted by switching the sign of inter- 
action from attractive to repulsive. The jets are distin- 
guished from the burst in that the directions of the jets 
arc purely radial. 

To understand this phenomenon, we performed nu- 
merical simulations under conditions similar to those of 
the experiment. Figure 7 shows the gray-scale images 
of the integrated column density after a t cxp and = 3.6 




radial direction 



Fig. 8. Snapshot of the collapsing BEC just before the collapse is 
interrupted. [Reproduced from Fig. 2 of H. Saito and M. Ueda, 
Phys. Rev. A 65, 033624 (2002).] 



ms expansion in (a) and 5.5 ms expansion in (c). These 
results look very similar to the experimental ones. 8 -* 

The parallel fringe pattern seen in Fig. 7 is reminis- 
cent of the interference fringe between waves emanating 
from two point sources. In fact, at the time when the 
interaction is switched from attractive to repulsive, the 
collapsing condensate exhibits two spikes in the atomic 
density that serve as two point sources of matter waves, 
as shown in Fig. 8. 

Thus the experimental observations of jets strongly 
suggests that atoms emanating from the spikes are co- 
herent. 

6. Pattern Formation 

We next consider a situation in which a BEC initially 
having a repulsive interaction is prepared and then the 
sign of interaction is suddenly switched to an attractive 
one. 

Figure 9 shows the time evolution of the integrated 
column density of the collapsing BEC. 5 ' 7 ) We initially 
prepare a repulsive BEC with fifty thousands atoms, and 
at time t = switch the sign of the interaction from 
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Fig. 9. Time evolutions of the integrated column densities of three 
collapsing BECs: (a) isotropic trap, (b) cigar-shaped trap, and 
(c) pancake-shaped trap. [Reproduced from Figs. 7 and 8 of H. 
Saito and M. Ueda, Phys. Rev. A 65, 033624 (2002).] 



repulsive to attractive. 

We can see that as time elapses, shells of high atomic 
density gradually form. At time t — 7.7ms, four shells 
are formed. But shortly after this formation (t = 8.2ms), 
the innermost shell collapses and disappears. Then, one 
by one, the other shells collapse and disappear as they 
reach the center of the trap, until finally all the shells 
have disappeared completely. The concentric circles in 
this gray-scale image clearly show the formation of the 
shell structure of the atomic density. This reveals a new 
self- focusing effect of the attractive BEC. 

For the cigar-shaped trap a cylindrical-shell structure 
is formed [Fig. 9(b)], while for the pancake-shaped trap 
a layered structure is formed [Fig. 9(c)]. 

7. Collapse of attractive BECs with a vortex 

Finally, we discuss the collapse of attractive BECs 
with a vortex. 17 ^ This subject is of interest because, gen- 
erally, attractive interactions cause the atoms to gather 



at the center, but in this case such gathering will be 
prohibited because the density right at the center must 
be zero due to the topological constraint of the vortex. 
We thus expect an interesting competition between the 
self-focusing effect and the topological constraint. 

Unfortunately, attractive BECs cannot hold vortices 
in any thermodynamically stable state. However, the 
Feshbach technique enables us to prepare a vortex in an 
attractive BEC by first creating a vortex in a repulsive 
BEC and then switching the sign of interaction to at- 
tractive. What are the ensuing dynamics of this system? 

To study this problem we work with the GP equa- 
tion (1) in two dimensions, and assume that the ini- 
tial state is a single vortex state. When the strength 
of attractive interactions exceeds a critical value, the 
Bogoliubov spectrum for the quadrupole mode acquires 
an imaginary part. This signals the onset of dynamical 
instability of the quadrupole mode. By increasing the 
strength of attractive interaction, we obtain a series of 
dynamical instabilities consisting of quadrupole, dipole, 
and monopole instabilities. 

Figures 10(a)-(f) show the time evolution of the den- 
sity and phase profiles of a vortex state with g = —9, 
which exceeds the critical value for the quadrupole mode, 
but does not exceed that for the dipole mode. Due to 
the quadrupole instability, the vortex splits into two clus- 
ters that rotate around the center of the trap, the two 
clusters reunite to recover the original shape, and this 
split-merge cycle repeats. 

The insets show the phase plots. At t = 16, three 
branch cuts appear. The central one exists from the 
outset, and the other two emerge as the vortex splits, 
in accordance with the fact that the m — 3 component 
grows upon the vortex split. The two side vortices can- 
not be seen in the density plot, and thus they may be 
considered "ghost vorticies" that carry a negligible frac- 
tion of the angular momentum. 

When the size of the system becomes too small to be 
observed by the in situ imaging, we may expand the con- 
densate by switching the sign of interaction to positive. 
Figures 11(g) and (h) show snapshots of such expanded 
images. We find that the stronger the repulsive interac- 
tions, the more fringes appear. 

Figure 11 shows the case in which g exceeds the criti- 
cal value for the dipole instability. The dipole instability 
gives rise to a transferring of the atoms from one cluster 
to the other, thereby inducing the collapse of the con- 
densate. In Figs, (a) to (c), cluster A grows, and then 
cluster B grows in a seesawing motion, until finally clus- 
ter B has absorbed most of the atoms and collapsed. In 
the case of (d)-(f), both clusters collapse immediately af- 
ter the vortex split. In this collapse process, we find the 
exchange of a vortex-antivortex pair. Figures (g) to (i) 
show the case with even stronger attractive interactions, 
in which the monopole instability appears. We see that 
the vortex first shrinks isotropically, then splits into two 
pieces, and finally collapses. 
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Fig. 10. Split-merge cycle of a vortex with attractive interaction 
(a)-(f) and expanded images (g) and (h). [Reproduced from Fig. 
3 of H. Saito and M. Ueda, Phys. Rev. Lett. 89, 190402 (2002).] 



8. Summary and Remarks 

We conclude this paper by summarizing our answers 
to each of the questions presented in the abstract. 

• BEC can survive partially after the collapse because 
the collapse occurs only in extremely localized re- 
gions and nowhere else in the BEC. 

• For an isotropic trap, the number of remnant BEC 
atoms is always below the critical number. How- 
ever, for an anisotropic trap the number exceeds the 
corresponding critical number, since after the inter- 
mittent implosions the shape of the gas becomes 
highly distorted and the collapse does not necessar- 
ily occur at the center of a BEC. 

• The atom burst is caused by a slight surplus of the 
zero-point kinetic energy over the attractive inter- 
action due to the removal of atoms in the highly 
localized high-density region by the three-body re- 
combination loss. Note that in the course of the 
collapse, the attractive interaction surpasses the ki- 
netic energy; the sudden removal of atoms due to 
the loss breaks the balance slightly in favor of the 
kinetic energy, leading to a mean burst energy of as 
low as lOOnK. 

• The gradual decay in the number of remnant BEC 
atoms, reported in Fig. 1 of Ref., 8 ) strongly suggests 
that the intermittent implosions occur during the 
decay process. 

• A sudden interruption of the collapse by switching 
the sign of the interaction from attractive to repul- 
sive leaves the collapsing BEC highly localized at 
two points for the cigar-shaped geometry under the 
conditions used in Ref. 8 ) These two localized high- 
density regions subsequently play the role of point 
sources of matter waves and yield the observed par- 
allel fringe patterns. This explains why the jets oc- 
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Fig. 11. Collapse of a vortex when the dominant unstable mode is 
the dipole mode (a)-(f) and the monopolc mode (g)-(i). [Repro- 
duced from Fig. 4 of H. Saito and M. Ueda, Phys. Rev. Lett. 
89, 190402 (2002).] 



cur only in the radial directions, unlike the bursts, 
and suggests that the jets are coherent. 

It is interesting to observe that the atom burst also 
occurs when the positive s-wave scattering length is sub- 
stantially increased in a very short time. 25,26 ) The dif- 
ference is that in these experiments the magnetic field is 
ramped up very close to the Feshbach resonance, where 
bound-state molecules are likely to be formed. In the 
Bosenova experiments, in contrast, the magnetic field is 
ramped away from the resonance, and therefore molecu- 
lar formation is expected to play a much smaller role than 
in the case of repulsive interaction. However, it seems 
likely that the two-body correlation is significantly en- 
hanced in the course of collapse. Moreover, in the midst 
of the collapse the peak density becomes so high that the 
quantum depletion of the condensate probably cannot be 
ignored. Nonetheless, our extensive study described in 
this paper suggests that a generalized GP equation (2) 
together with the formula (3) offers a consistent overall 
picture of the Bosenova experiments. 
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